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£N| , Abstract. The NA60 experiment measured dimuon production in In-In collisions at 

158 AGeV. This paper presents a high statistics measurement of <f> — ► ji^i with the 
^ . specific objective to provide insight on the (f> puzzle, i.e. the difference in the inverse 

1} \ T slopes and absolute yields measured by NA49 and NA50 in the kaon and lepton 

O i' channel, respectively. Transverse momentum distributions were studied as a function 

of centrality. The slope parameter T shows a rapid increase with centrality, followed 
by a saturation. Variations of T with the fit range of the order of 15 MeV were 
observerd, possibly as a consequence of radial flow. The 4> meson yield normalized to 
the number of participants increases with centrality and is consistently higher than the 
yield measured by the NA49 experiment at any centrality. 

(N : 
^t- ■ 

One of the most striking predictions of QCD is the occurrence of a phase transition 
from hadronic matter to a deconfmed plasma of quarks and gluons, when sufficiently 
high energy densities are reached. Strangeness enhancement was proposed long ago as 
a signature of such phase transition [1] . Experimentally, a rather interesting observable 
is the meson, due to its ss valence quark content. Moreover, it has been argued that 
close to the phase boundary, the spectral function of the 0, and thus its mass, width and 
branching ratios, could be modified [2]. The yield and transverse momentum spectra 
of the meson were studied extensively at the SPS in Pb-Pb collisions at 158 A GeV. 
The NA49 experiment measured the — > KK channel [3] (with statistics limited to 
Pt < 1-6 GeV), while the NA50 experiment measured the — > fifi channel [I] (with 
acceptance limited to pt > 1-1 GeV). Both experiments observed an enhancement of 
the yield with the size of the collision system. However, the absolute yields measured 
in the common kinematic window (transverse momentum, rapidity and Collins-Soper 
angle) disagree by a factor of about two [5]. It was suggested that kaon absorption 
or rescattering prevents the reconstruction of most in-matter — > KK decays, in 
particular at low transverse momentum, while the mesons decaying in the lepton 
channel would not be affected [21 El 13- The effect predicted by these models, however, 
is smaller than the observed difference. NA49 and NA50 disagree also on the inverse 
T slopes determined from the transverse momentum spectra. The NA49 inverse slopes 
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Figure 1. T slope parameter as a function N par t compared to the results of NA49 
and NA50. Left: fit performed in the NA49 range (0 < p T < 1.6 GcV). Right: 
Fit performed in the NA50 range (1.1 < pt < 2.6 GeV; due to low statistics, the 2 
most peripheral points are integrated). Square brackets: total error (statistical and 
systematic). 



are consistently larger and with a stronger centrality dependence than those measured 
by NA50. This could be partly due to radial flow which tends to make the T slopes 
extracted from a low pt region higher. A further increase of the NA49 slopes might be 
due to the depletion of the — > KK at low transverse momentum. The discrepancy 
in absolute yields and slopes became known as the puzzle [2J [8] . Production of in 
central Pb-Pb collisions was recently studied also by the CERES collaboration [9] both 
in the kaon and in the dielectron channel. The measured inverse slope and yield are 
compatible with the NA49 results in both channels, although the ee measurement is 
affected by large uncertainties. 

This paper presents a high statistics measurement of — > fifi in In-ln collisions 
at 158 A GeV, with the specific objective to provide insight into the puzzle. The 
sample consists of 360000 signal pairs, with ~70 000 events lying in the peak. The 
signal/background ratio integrated over centrality is ~ 1/3 below the peak. The data 
were divided in 5 centrality bins, selected on the basis of the number of participants, 
with (Npart) = 15, 39, 75, 132, 183. The number of participants was estimated from the 
charged tracks multiplicity measured in the vertex detector, matching the knee at high 
multiplicities to a Glauber calculation and assuming a linear relation N part cx N c h- For 
more details on the analysis procedure and on backgrounds subtraction see Ref. [TOlfTT]. 

In order to study the px distributions, events in a small mass window centered at 
the pole were selected (0.98 GeV < M w < 1.06 GeV). The spectrum of the continuum 
below the was estimated and subtracted selecting events in two side mass windows 
(0.88 GeV < M wt < 0.92 GeV and 1.12 GeV < M w < 1.16 GeV). After subtraction 
of the combinatorial background, physical continuum and fake tracks, the kinematic 
distributions were corrected for acceptance times reconstruction efficiency, estimated 
with an overlay Monte Carlo simulation (a Monte Carlo dimuon is reconstructed on 
top of a real event). In principle this correction should be evaluated as a 3-dimensional 
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matrix, as a function of pr, y and cos 6cs- This approach also requires fiducial cuts 
to remove the phase space corners, with a corresponding loss in statistics. In order to 
preserve the available statistics, a 1-dimensional correction as a function of the variable 
under study was used, after a careful tuning of the input kinematic distribution via 
an iterative procedure. In order to extract the inverse slope parameter T, the pr 
distributions were fitted with the exponential function l/p T dN/dpT = exp (— m^/T). 
The NA60 acceptance goes down to = GeV and the pt spectra have statistics up to 
> 2.6 GeV. The fit was repeated in the full range 0.0 GeV < pt < 2.6 GeV and in the 
two sub-ranges (0.0 GeV < < 1-6 GeV and 1.1 GeV < < 2.6 GeV), corresponding 
to the experimental windows of the NA49 and NA50 experiments, respectively. The 
X 2 /ndf of the fits ranges between 0.6 and 1.5. The T measured in the full range shows a 
rapid initial increase with centrality, followed by a saturation. It should be noticed that 
in presence of radial flow the T slopes extracted from exponential fits depend on the fit 
range [12]. The results in the sub-ranges are shown in figure [fl compared to the results of 
the other experiments. As seen in the figure, when the fit is performed at high p?, there 
is a decrease of the average value of T and a flatter trend with centrality. The maximum 
difference, in the most central bin, is of the order of ~ 15 MeV. The discrepancy between 
NA50 and NA49 seems not consistent with the T slope variation observed in the NA60 
data, suggesting that some other mechanism like kaon absorption/rescattering may be 
at play. A direct comparison of the absolute values of T between NA60 and the Pb-Pb 
experiments is not straightforward due to possible system-size effects. 

Two different approaches with largely independent systematics were used to 
determine the average number of mesons produced per interaction, (0). In the first, 
the cross section integrated in centrality has been evaluated using the observed 
number of mesons and the number of incident ions measured with 2 independent beam 
counters, ((f)) was then calculated dividing by the total inelastic In-In cross section, 
estimated with the general formula Omei^AB = nrl [A 1 / 3 + B 1 / 3 — j3 (A^ 1 / 3 + B^ 1 / 3 )] , 
where A and B are the mass numbers of the two nuclei, r = 1.25 fm and (3 = 0.83 [13J. 
In the second approach the J/ip was used as a reference. (0) can be calculated from the 
measured ratio (0) / (J/ip) = (N^/Acc^ x e rec ^) / (Nj/^/ Accj/^ x e reCi j/^). Once corrected 
for anomalous and nuclear absorption, (J/ip) scales with the number of binary collisions 
and can be written as (J/ip) = (o~nn/o~nn) • N co u. By multiplying this value by the 
above ratio, one obtains (0). The 2 methods agree within ~ 10%. The enhancement 
is determined by normalizing the yield to the number of participants. Figure [2] shows 
the ratio (0) /N part and the comparison to the NA49 measurements in several collision 
systems [T4"[ [TBI [16] . The enhancement in In-In with respect to p-p reaches a factor 
of five in the most central collisions. As seen, the observed enhancement observed by 
NA49 is consistently smaller at any centrality. NA49 showed that the — > KK breaks 
Np ar t scaling, and that (0) /N part for central collisions and for several collision system 
already saturates for Si-Si collisions. The value of (0) /N part for central In-In collisions 
in the muon channel exceeds the corresponding values for Pb-Pb collisions in the kaon 
channel, suggesting that the NA60 \x\i data do not follow the trend measured by NA49 
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Figure 2. <p yield normalized to the number of participants in full phase space and 
corrected for the branching ratio. Compared to the NA49 measurement in several 
collision systems. The error bar of the Pb-Pb data includes statistical and systematic 
error. Square brackets: total error (statistical and systematic). C, Si and peripheral 
In points displaced by ±2 participants to improve readability. 



in the kaon channel. A further reduction of the systematic uncertainties is currently 
in progress. An unambiguous comparison to the NA50 result in full phase space is 
not possible, due to lack of consensus on the value of the slope parameter in Pb-Pb 
collisions (see above). An extrapolation with the extreme hypotheses T = 220 MeV and 
T = 300 MeV leads to much larger values with respect to In- In, even in the upper-bound 
scenario T = 300 MeV. The measurements of the CERES experiment do not rule out 
differences of the order of those observed by NA60, due to the large uncertainties on the 
ee measurements. 

In conclusion, the long-standing <p puzzle is yet to be clarified. The on-going analysis 
of the — > KK channel in the NA60 data may be decisive to clarify the picture. 
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